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181. 


ON SIR W. R. HAMILTON’S METHOD FOR THE PROBLEM OF 
THREE OR MORE BODIES. 


[From the Quarterly Mathematical Journal, vol. 11. (1858), pp. 66—73.] 


THE problem of three or more bodies is considered by Sir W. R. Hamilton in 
his two well-known memoirs on a general method in Dynamics, Phil. Trans. 1834 
and 1835, and the differential equations for the relative motion, with respect to the 
central body, of all the other bodies are obtained in a form containing a single 
disturbing function only. Several methods of integration are given or indicated, among 
others, one which is in fact the method of the variation of the elements as applied 
to the particular form of the equations of motion. But the investigation shows (and 
Sir W. R. Hamilton notices this as a defect in his theory, as compared with the 
ordinary theory of the variation of the elements), that in the method in question, the 
elements are not osculating elements, ie. that the positions only, and not the velocities 
of the bodies, can be calculated as if the elements remained constant during an 
element of time. The peculiar advantage of the method is of course the having a 
single disturbing function only, and this seems so important, that. if I may venture to 
express an opinion, I cannot but think that the method will ultimately be employed 
for the purposes of Physical Astronomy. But, however this may be, it has appeared 
to me that it may be useful to present the method in a separate and distinct form, 
disengaged from the general theory as an illustration of which it was given by the 
author; and this is what I propose now to do. 


Consider a central body M, and two other bodies M,, M,, and let the coordinates 
of M referred to a fixed origin be a, y, z, and the coordinates of M,, M, referred to 
the body M as origin be æ, Jı, 4% and a, Yz, Za respectively. Then the coordinates 
of M,, M, referred to the fixed origin, are +m, y+, 2+42, and æ +2, y+42, 2+2, 

iO: TH: 13 


www.rcin.org.pl 


98 ON SIR W. R. HAMILTON'S METHOD FOR THE [181 


respectively, and if as usual 7 denotes the Vis-viva or half sum of each mass into 
the square of its velocity, and U denote the force function, then we have 
T= 4M (#2*+y%4+2%), 
+ 4M, (£ +Y HU HYH aY), 
+ 4M, {(a! + af +y + yo’ + (2 + zY}, 
MM, 
y (a? + Yi + 2?) 
MM, 
y (x? + Yo a z2) 


PEE ES a a a 
N (a p La)" + (Y = Yo) + (2 = 23)? i 


Uz 


-|- 


and the equations of motion are as usual 


dd? _av_av 
dt dx’ dx dæ’ 
&e. 


If we assume that the centre of gravity of the bodies is at rest, then we have 
Ma’ + M, (x + a’) + M, (æ + a/)=0, &e., 


and consequently 


of = — ha + Maey r & My + My J = Mz’ + Miz! 
M+M, +M, I M + M, + M,’ M+M, +M,” 


Now the value of T is 
T= (M+M, + M)(a + y+ 2%) 
+ 2! (Mya,' + Mær) + y (My + Mey’) + 2 (Mizi + Maz) 
+ $M, (a? + yp’ + 4") 
+ $M, (a7 + y: + 22"), 
or, putting for a’, y’, 2 their values, 
T= $M, (@?+y2+2") 
+ $M, (21? + ys? + 23°) 


1 r t / , 1 
a HTE FE. (Mæ; + Mæ} + (My + Maye’? + Ma! + M:2Y}, 
and with this new value of T the equations of motion still are 


dt dx da, da,’ &e. 
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Suppose now that the differential coefficients of T, with respect to a’, Y1, 2; a, Yo, Za, 
are respectively P,, Q, Ri; Pa Qo Re, ie. write 


and imagine T expressed as a function of P,, Q, Ri; Pa Qa R, and when this is 
done put H=T—U (so that H stands for a function of P,, Q, R; Pa, A POR A 
Li, Yis 2; La Yas 22), then the equations of motion in Sir W. R. Hamilton’s form are 


ae E ai 
Now from the last given value of T 
Þi=Ma -pA (Mhe + Mn) 
1 = 44,2, M+ M, +M, 1% 2% ), 
P, = S A — op — — (Ma! + Me) 
2 2⁄2 M + M, + M, 1 22 /) 
and thence 
M ' ; 
Popp EI EI Ah + Mær), 
and consequen tly 
Mæ = Pia T (P, +P), 
, M, 
Mn, = P, + -q (Pit P) 


and we have 


- ie $ T P+ PY T {0 He + af +1R, + = (R, + R | 


PA 
2M, 


+ guy Pr Her Pah e H+ Oo} {m+ gr ces 20} 


= stg (M+ M, + M3) (P + P,P + (Qs + Q) + (B+ BM, 


or, reducing, 


T= (stay) Prt Qr+ Ro 
+ (amr, + am) PP + OF + BO 


1 
+ y (PiPs+ GQ. + RR), 
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and consequently x 
H= (54+ aq) (P+ O+ Re) 


+ (srt am) (Pè + Q+ Re) 


1 
+ wih: T QQ: + RR) 


y (a2 + Y? ig 2°) 
y (a2 + Ys + Z 2) 
= MM, 
y {(a, — t) + (% = Yo) +(4,- 22)*} , 
aud H having this value, the equations of motion are as before mentioned 


do, dH dP, dH , 
Per? A dat 


Instead of H write H+T where 


1 1 5 
1 i 
MM, 
(a? + yy rs 2?) 
MM, 
y (æ? + yf + 2°) i 


and 
T= y (PiPa+QQi+ RR) 
fe M,M, 
V{(@ — t} + (Yr — Yo)? + (A — 2}? 


and the function T is to be treated as a disturbing function. The equations of 
motion for the body M, become 


de, M+Mp at dP, ___ MMe, _dT 
dt MM, ` a q) idi- (a@P@tyeta)t day’ 
dy, _ fs ce dQ, MM dT 
Q+ a a i FOr ge Oe 
dt z dt (metye +z dh 
w nee he: Men cd ee 
dt MM, dR,’ dt (£e +y? +2 dz,’ 


and there is of course a precisely similar system of equations of motion for the body M,. 
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If we neglect T, the left-hand equations show that P,, Q, R, denote the 
velocities or differential coefficients sd ; Sir os multiplied by the constant factor 
re , and substituting these values in the right-hand equations, we obtain the 

1 


ordinary equations for the elliptic motion of the body M,; and similarly for the body 
M, We may, if we please, complete the solution by the method of the variation of 
the arbitrary constants. Suppose for this purpose that a, b, a, 4, fi, g9 are the 
elements for the elliptic motion of the body M,, then treating these elements as 
variable we must have 


dr, da, dx, db, , de, dp av 
da, dt + ab, dt coe t dg, dt = dP,’ &e., 


dP, da, , dP, db, dP, dy, dT 


da, dt db, dt” dg, dt ~~ day’ 

and it appears from these equations that as already noticed the disturbed values of 
the velocities are not (as they are in the ordinary theory) identical with the undisturbed 
values. 


The disturbing function T may be considered as a function of the elements of 
the two orbits and of the time, and it is easy to obtain, as in the ordinary theory, 
da 


the values of the differential coefficients Te’ &c. in the form 


da _ dT dT dT 
de = (™ bs) 7, + (a, Oaa oe #(G, A) Ty, > 
where 


Ela, b) , èla Bb) , 8(a, b) 
(a; ete FPS PY Ò (J, Qi Be, B)’ 


5(a,, b,) da db da, db, 


ò (@, Pa) da; dP, dP, da: 


if for shortness 


It will be remembered that in the ordinary theory, if 2 denote Lagrange’s dis- 
turbing function (Q =— R if R is the disturbing function of the Mécanique Céleste) 
the corresponding formulæ are 


(a b) + (a, oO .. +@ DTe 
where 
8 (a, b) 
Ò (£, x) 


S(a, b)  8(a, b) 


Go Sy, y) 5,2)’ 


= 


if for shortness 
S(a, b) da db da db 


S(@, 2) da’ de da da’ 
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or, what is the same thing, where 
_ 8, b) (a, b) _ 8a, c) © 
ola a) “Oty, ¥) (2, #)’ 


8(a, b) _da db _ da db 
Ò (æ, £) dada’ da’ da’ 


(a, b)= 


and 


Now the values of the coefficients (œ, bı), &. depend merely on the form of the 
expressions for a, b, &. in terms of P,, Q, R, m, Yı» 4 and t; hence comparing 
the two systems of formule and observing P,, Q, Æ (which in the formule for the 
present theory correspond with æ, y, 4 in the other system of formulz) are 
respectively equal to æ’, Yy, 2%, each of them multiplied by the constant factor 
MM, 
M+M,’ 
of elements in the present theory are at once deduced from- the formulæ for the 
variations of the same system of elements in the ordinary theory by writing —T in 
_the place of Q and multiplying the values of the variations by the constant factor 
MM, 

M+M,’ 


Take then as elements Jacobi’s canonical system (?), viz. if we put 


it is easy to see that the formule for the variations of any given system 


a, the semiaxis major, 

@, the eccentricity, 

@,, the longitude in orbit of pericentre, 
&, the mean longitude in orbit at epoch; 
6,, the longitude of node, 


¢:, the inclination, 


m, the mean motion f- vi (= > =! : 


then the canonical elements are 


and 


A = -— fma’, 
$, =  maN(l-e’), 
G = man(l-—e’) cosg, 


BA = + (G = @), 
©, = = A, 
H = 0, > 


1 I have for uniformity adopted Jacobi’s canonical system, see his paper ‘‘Neues Theorem der analytischen 
Mechanik,” Crelle, t. xxx. pp. 117—120 (1846); but it is proper to remark that Sir W. R. Hamilton, in his 
Memoirs above referred to, employs a slightly different but equally elegant system of canonical elements, 
and that the discovery of such a system belongs to Sir W. R. Hamilton, and is part of his general theory. 
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(the signs of the two elements A, §, have been changed, but this makes no difference 
in the formule) then the equations for the variations of the elements are 


dA M + M, dT 


dt MM, d’ 
qB,__M+M, dt 
al MM, d®&,’ 
d@,_ M+M, aT 
a GMM, dS,’ 


q,_ , M+ M, aT 
dE = MM tg 


d®, M+M, dT 


oe. me. oe 
dt * MM, dG,’ 


and it is easy thence to deduce the formule for the variations of any system 
of elements which it may be thought proper to make use of, for instance the system 


a, ĉi, Bi, €, 0, Qı. 


It will be recollected that in the preceding system of formulæ the value of the 
disturbing function T is 


T= 7 PsP: + QQ + RiR) 


Hwy ops) atah Tae oomi: 
V{(@, — 2a} + (1 — YF + (a a) 
and that as a first approximation P,, Q,, R, are respectively equal to the velocities 


zy, h, Z, each multiplied by the constant factor weir and P,, Q, R, are 
1 


respectively equal to the velocities æy, Y7, 2%, each multiplied by weir , 
2 


2, Stone Buildings, 18th Oct., 1856. 
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